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ABSTRACT 



We analyze the dynamics of gas-dust coupling in the presence of stellar radiation pressure in 
circumstellar disks, which are in a transitional stage between the gas-dominated, optically thick, 
primordial nebulae, and the dust-dominated, optically thin Vega- type disks. Meteoroids and dust 
undergo radial migration, either leaving the disk due to a strong radiation pressure, or seeking 
^ ' a stable equilibrium orbit in corotation with gas. In our models of A-type stars surrounded by 

a total gas mass from a fraction to dozens of Earth masses, the outward migration speed of 
dust is comparable with the gas sound speed. Equilibrium orbits are circular, with exception of 
those significantly affected by radiation pressure, which can be strongly elliptic with apocenters 
extending beyond the bulk of the gas disk. 

The migration of dust gives rise to radial fractionation of dust and creates a variety of pos- 
sible observed disk morphologies, which we compute by considering the equilibrium between the 
dust production (at radii < 10AU) and the dust-dust collisions removing particles from their 
O equilibrium orbits. Sand-sized and larger grains are distributed throughout most of the disk, 

with concentration near the gas pressure maximum in the inner disk. Smaller grains (typically 
■ in the range of 10 to 200/im) concentrate in a prominent ring structure in the outer region of 

the gas disk (presumably at radius ~100 AU), where gas density is rapidly declining with radius. 
The width and density, as well as density contrast of the dust ring with respect to the inner dust 
disk depend on the distribution of gas and the mechanical strength of the particles, but do not 
depend on the overall dust production rate. 

Our results open the prospect for deducing the distribution of gas in circumstellar disks by 
observing their dust. We have qualitatively compared our models with two observed transitional 
disks around HR 4796A and HD 141569A. (Gas component has been detected, but not yet 
mapped in detail, in the second object.) Dust migration can result in observation of a ring or 
a bimodal radial dust distribution, possibly very similar to the ones produced by gap-opening 
planet(s) embedded in the disk, or shepherding it from inside or outside. We conclude that a 
convincing planet detection via dust imaging should include specific non-axisymmetric structure 
(spiral waves, streamers, resonant arcs) following from the dynamical simulations of perturbed 
disks. 

Subject headings: accretion, accretion disks — circumstellar matter — planetary systems: formation 

1. Introduction quence stars. Because their evolution is intimately 

connected with the process of planet formation, 
Dust and gas disks surround stars of different studies of dust amountj size spe ctrum, mineral- 

ages, from pre-main sequence to post-main se- ogyj and spatial distr ibution in disks can provide 
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valuable information about the crucial stages of 
planet formation. 



1 



The amount of dust in disks generally decreases 
with time (Zuckerman & Becklin 1993; Natta, 
Grinin, & Mannings 2000). As a consequence, 
the fraction of detectable, optically thick disks 
(at <3 AU from the star), decreases from nearly 
100% at the estimated object age ~0.3 Myr down 
to a few percent at the age of ~10 Myr (Hillen- 
brand & Meyer 1999; Meyer & Beckwith 2000). 
The amount of gas decreases as well (Zuckerman, 
Foreville, & Kastner 1995; Liseau & Artymow- 
icz 1998), especially abruptly at an age of sev- 
eral to ~10 Myr for solar-type stars, when cir- 
cumstcllar disks lose most of the primordial proto- 
stcllar/protoplanetary material through processes 
such as viscous accretion and photoionization, as 
discussed by Hollcnbach et al. (2000). From the 
original value of close to 100, the gas-to-solids (if 
not necessarily gas-to-dust) mass ratio in the T 
Tauri and Herbig Ae/Be stars decreases during 
this period (i.e., typically at age 3-10 Myr) to 
much smaller values, eventually much less than 
unity, thus leaving a predominantly dusty remnant 
disk. 

Vega-type stars, observable owing to their large 
thermal infrared (IR) excess over the extrapolated 
IR stellar photospheric radiation, and sometimes 
resolved in scattered visible/near-IR starlight, are 
thought to contain such gas-poor disks. The prime 
example of an object of this kind is (i Pic, which 
is 20-100 Myr old, and has a large disk (>10 3 AU 
radius) containing gas and dust, apparently in a 
ratio smaller than unity (Artymowicz 1997; Liseau 
& Artymowicz 1998). Based on dust destruction 
timescale in Vega- type disks, the disks are thought 
(or sometimes defined) to be continuously replen- 
ished by the collisions and evaporation of plan- 
etesimals. Most convincingly in (3 Pic, but also 
in a few younger, spectroscopically variable Her- 
big Ae stars, the existence of planetesimals and 
planetary objects necessary to perturb planetes- 
imal orbits, has been inferred from spectroscopy 
(Lagrange, Backman, & Artymowicz 2000; Grady 
et al. 2000). 

The disks around stars with ages between 3 
and 10 Myr are particularly interesting, because 
of the expected presence of protoplanets, grow- 
ing by collisional accumulation of planetesimals 
and through gas accretion (Hayashi, Nakazawa, 
& Nakagawa 1985; Lissauer 1993). We shall call 
such disks transitional, in reference to a transition 



from a relatively high-mass, optically thick pri- 
mordial nebulae to low-mass, optically thin dust 
disks. (Other names, related to the central star, 
are used for the objects in which such disks re- 
side: young main sequence stars, or old pre-main 
sequence stars; cf. Lagrange et al. 2000.) 

Recent imaging of two apparently transitional 
disks, HR 4796A and HD 141569A, provided in- 
teresting insight into the spatial distribution of 
dust. The AO V star HR 4796 A with age 8 ± 2 
Myr (Stauffcr, Hartmann, & Barrado y Navas- 
cues 1995) has been imaged, first in mid-IR and 
afterwards in near-IR (Jayawardhana 1998; Ko- 
erner et al. 1998; Schneider et al. 1999; Telesco 
et al. 2000). The second object, B9.5 Ve star 
HD 141569A, is 2-10 Myr old. Near-IR images 
of its dust disk were obtained by Weinberger et 
al. (1999), and Augereau et al. (1999a), while 
the thermal mid-IR flux was mapped by Fisher et 
al. (2000). The disk in the image of HR 4796A 
extends out to about 100 AU, a distance 4 times 
smaller than in HD 141569A, and is significantly 
narrower than its counterpart (of order 20 AU, vs. 
200 AU). The two systems have, however, very 
similar dust covering factor, or the ratio of to- 
tal IR dust luminosity to the stellar luminosity, 
/dust = L dust /L„, equal to 0.7% in HR 4796A, 
and 1% in HD 141569A, which are a factor of 3 to 
4 times larger than in (3 Pic. 

Both disks show structure, which can be mod- 
eled as axisymmetric, broad rings with inner holes 
(near side-far side asymmetry of surface brightness 
by a factor of 1.5 or less, has been interpreted as 
being due to an oblique viewing angle combined 
with anisotropic scattering by dust.) The radii 
of the inner holes are ~ 60 AU for HR 4796A 
(Kocrncr et al. 1998; Schneider et al. 1999) and 
- 150 - 200 AU for HD 141569A (Augereau et 
al. 1999a; Weinberger et al. 1999). The image 
of HD 141569A taken at 1.1/im by Weinberger et 
al. (1999) has been analyzed to reveal a double 
ring, or ring-gap-ring morphology with the shal- 
low gap (or dip) in the surface brightness, cen- 
tered on the projected distance r = 250 AU. The 
two maxima are found at 200 and 350 AU. Note- 
worthy, the same system observed at 1.6/im by 
Augereau et al. (1999a) shows only one broad 
ring peaking in brightness at r = 325 AU. It is 
not clear to what extent the apparent differences 
in the inferred surface brightness are due to the pe- 
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culiarities of observational techniques, such as the 
different time lag between the registration of ob- 
ject and comparison star (point-spread function), 
and to what extent by intrinsic factors like, e.g., 
the possible wavelength-dependent scattering effi- 
ciency of dust around HD 141569A. For example, 
if the inner peak of HD 141569A is dominated by 
small grains, then it would Rayleigh-scatter 4.5 
times better at l.l^m than at 1.6/xm. Recently, 
HD 163296, a Herbig Ae star, was announced by 
Grady et al. (1999) to have a disk-and-gap mor- 
phology similar to that of HD 141569A. 

Disks around Vega- type stars, which are more 
evolved than the transitional disks discussed 
above, also have inner holes or clearing regions, 
mostly but not always completely devoid of dust 
(for review see Lagrange et al. 2000). Mid-IR im- 
ages of the (3 Pic disk show the central low-density 
hole (Lagage & Pantin 1994; Pantin, Lagage, & 
Artymowicz 1997) and submillimeter images of 
Fomalhaut and e Eridani show the ring-like struc- 
ture of dust disks (Holland et al. 1998; Greaves 
et al. 1998). It is expected that the gravity of 
embedded large bodies, such as planets or brown 
dwarfs, form such inner holes in the dust distribu- 
tions (Roqucs et al. 1994: Liou & Zook 1999). 

The inner holes of the transitional disks around 
HR 4796A and HD 141569A may also be produced 
by planetary bodies embedded in the disks (or 
rather, in these holes). Moreover, planetary influ- 
ences creating the disk "gap" at 325 AU from HD 
141569A (and a similar distance from HD 163296) 
have been proposed as an explanation of these 
striking features. Weinberger et al. (1999) es- 
timate that a 1.3 Jupiter mass perturber would 
open a gap with the observed extent. Grady et 
al. (1999) conclude that their observations lend 
strong support for early formation of massive plan- 
ets at distances as large as r > 300 AU from the 
stars. 

However, there arc several potential difficulties 
with the proposal that giant planets form rou- 
tinely at such large distances in disks, in only a few 
Myr. Notice that the circular orbit of a planet, re- 
quired by the circularity of a narrow "gap" , would 
imply a local formation history of the hypothetical 
planet. Formation in the inner disk and ejection 
to a large distance by, e.g., planet-planet pertur- 
bations would not provide a nearly circular orbit 
required. Admittedly, radial migration of planets 



via interaction with the disk may also occur but 
an outward migration distance is always moder- 
ate. Standard theory of planet formation requires 
a high surface density of solids for giant planet for- 
mation at large radii. For example, the numerical 
simulations of planetesimal accretion at 70 AU by 
Kenyon et al. (1999) showed that the Pluto-sized 
planets (but not the required multi-Earth mass 
planets) may form in 10 Myr only if the mass of 
the disk is as large as 10 — 20 times the minimum 
solar nebula mass (M g <~ 0.2 — O.4M ), a mass 
about 10 times larger than the typical masses of 
disks around T Tauri stars (Osterloh & Beckwith 
1995). 

The amount of gas on the outskirts of a typical 
T Tauri disk, at the age of ~1 Myr, is probably 
also insufficient for other modes of giant planet for- 
mation, such as the disk fragmentation, although 
such a possibility in principle exists. One argu- 
ment against it, however, is apparent in the HD 
141569A system, where the surface brightness pro- 
file does not show any deep depletion of dust (an 
empty gap) such as those which quickly appear 
around a giant planet in hydrodynamical models 
(e.g., Lubow, Seibert, & Artymowicz 1999). On 
the contrary, the "gap" is merely a 14% dip in the 
surface density of dust, as compared with neigh- 
boring disk regions. Weinberger et al. (1999) thus 
propose a very recent timing for the formation of 
their proposed Jovian perturber. Troublesome as 
it is for any formation scenario, this special timing 
requirement seems particularly incompatible with 
the early disk fragmentation scenarios. 

Motivated by the apparent axisymmetry or 
near-axisymmetry of the observed ring and gap 
structures in disks, we devote this paper to the 
exploratory study of an alternative mechanism 
for the appearance of a variety of possible dust 
morphologies including outer dust rings and ap- 
parent disk gaps, connected with the phenomenon 
of radial migration of dust in gas disks. In the 
context of transitional disks, which are optically 
thin to the stellar visible and near-UV radiation 
and have gas component with moderate mass, we 
consider two important effects: radiation pres- 
sure from the central star and gas drag. As we 
will show, the combination of radiation and drag 
forces results in an outward direction of migra- 
tion for the observable dust and strongly affects 
the observed disk morphology, while preserving 
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the usual inward direction of migration of sand 
and pebbles. It is easy to see why transitional 
disks support the most vigorous migration. Pri- 
mordial protoplanetary nebulae are both optically 
thick and so dense (gas:dust mass ratio ~ 10 2 ) 
that most dust grains "freeze" in gas, and are not 
able to migrate through the gas. Vega- type stars, 
on the other hand, are dust dominated (Lagrange 
et al. 2000). They generally contain so little gas 
(gas:dust <C 1) that the motion of dust is weakly 
affected in the grain's lifetime. In contrast, HR 
4796A and HD 141569A may have gas disks com- 
parable with or somewhat more massive than the 
dust disks. Zuckerman et al. (1995) have found 
CO line emission from the gas inside 130 AU dis- 
tance from HD 141569A. A double Ha line from 
the gas disk in the immediate vicinity of the star 
(Dunkin, Barlow, & Ryan 1997) also demonstrates 
the existence of an extended gas disk, required for 
sustained accretion of gas. Whether this gas is 
a remnant of primordial gaseous nebula or a sec- 
ondary gas released by solids (for instance, evap- 
orating planetcsimals), is not known at present. 
In any case, transitional disks may have substan- 
tial amounts of gas (many Earth masses), even if 
an unfavorable viewing geometry of most systems 
precludes a spectroscopic detection of the circum- 
stellar absorption lines. Artymowicz (2000) argues 
that gas must be present and gas drag must be 
operative in the transitional disks, or else the ra- 
diation pressure-boosted destruction rate of solids 
would be prohibitively large. 

One feature of our model, which may be observ- 
able, is the radial fractionation of dust according 
to grain size. The most important feature is that 
at the outer edge of the gas disk, dust forms a 
ring without the help of any shepherding bodies. 
Therefore, the ring structure of the dust disk does 
not necessarily indicate the existence of planets or 
brown dwarfs. Klahr & Lin (2000) also proposed 
a model for the ring formation due to the dust mi- 
gration. Whether or not these models can quanti- 
tatively reproduce all the observational data is, of 
course, a matter of detailed modeling, taking into 
account each system's scale, central star's prop- 
erties, indications of dust mineralogy from spec- 
troscopy etc. We defer such detailed modeling of 
transitional disks to a later publication, but pro- 
vide qualitative discussions in §7. It should, how- 
ever, be stressed that our model is testable against 



the existing observations, and can suggest new ob- 
servations. Excluding dust migration as a cause of 
the observed disk structure might, in fact, be as 
important as finding that it matches the obser- 
vations of a given object, since such an exclusion 
would very strongly support the hypothesis of dy- 
namical sculpting by planets. 

The plan of the paper is as follows. In §2, 
we describe the equations of the dust grain mo- 
tion. In §3, we derive approximate expressions of 
radial velocity of grains. We show both inward 
and outward migrations of grains occur due to the 
combined effect of radiation pressure and gas drag 
forces. In §4, we show that grains segregate in the 
gas disk according to their sizes and that grains 
of medium size (10 — 100/im in our models) con- 
centrate at the edge of the gas disk. In §5, the 
orbital evolution of grains is calculated numeri- 
cally. We discuss the excitation of eccentricities 
of grain orbits at the edge of the gas disk. In §6, 
we calculate the lifetime of grains by taking into 
account the collisional destruction of grains. We 
derive the density structure of the dust disk and 
discuss the formation of an outer dust ring. In 
§7, we show a sample image in near- and mid-IR 
from our simple model. Finally, §8 discusses the 
relevance of the migration-induced dust structures 
to observed systems, and compare them with dust 
features of other origin. 

2. Equations of Motion of Dust Grains in 
a Gas Disk 

We consider the motion of dust grains in a gas 
disk. The equation of motion of a dust grain with 
mass m and the position vector r is 



GMm 



f + F rad + F PR + F g 



(1) 



where r is the distance of the grain from the cen- 
tral star, r is the unit vector directed to the grain, 
G is the gravitational constant, M is the cen- 
tral star's mass, and F rac i, -Fpr, and F g are the 
radiation pressure, Poynting-Robertson drag (PR 
drag), and gas drag force, respectively. 

The radiation pressure force F ra( j is expressed 
using the ratio of radiation pressure force to the 
gravitational force /?, 



n GMm^ 
Frad = P 



(2) 
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The ratio (3 is written as (Burns, Lamy, & Soter 
1979) 

= ZLQpr n) 
H WirGMcspd ' 1 ' 

where L is the central star's luminosity, c is the 
speed of light, s is the grain's radius, pd is the 
material density of the grain, and Qpb, is the radi- 
ation pressure coefficient averaged over the stel- 
lar spectrum. The values of (3 are calculated 
using Mie theory (Artymowicz 1988) and shown 
in Figure 1. The grains are assumed to con- 
sist of silicates. We use the laboratory measure- 
ments of optical constants of olivine (MgFeSiCU) 
and magnesium- rich pyroxene (Mgo.sFeo^SiOa) 
by Dorschner et al. (1995) 2 . The grains are mod- 
eled as porous, containing 50% volume fraction 
of vacuum (the volume fractions of the two sili- 
cates being 25% each). The density of composite 
grains is then p d = 1.25 g cm~ 3 . For the central 
stars, we used parameters of HR 4796A (solid line) 
and HD 141569A (dashed line), and the renormal- 
ized spectrum of AO V type star Vega. The mass 
and luminosity of HR 4796 A are M = 2.5M and 
L = 21.0L Q , respectively (Jura et al. 1993; Jura 
ct al. 1998), while for HD 141569A these param- 
eters are M = 2.3M Q and L = 22.4i Q (van den 
Anckcr ct al. 1998). 

The PR drag force Fpr is 



PR 



= -/3 



GMm 



V c c/ 



(4) 



where v is the velocity vector of the grain, and v r 
is the velocity component in the r direction (Burns 
et al. 1979). 

If the velocities of the grain (v) and the gas (v g ) 
differ, the grain experiences a gas drag force F g . 
The sizes of grains considered in this paper arc 
much smaller than the mean free path of the gas 
molecules. Thus, if the velocity difference A.v = 
v — v g is much smaller than the mean thermal 
velocity of gas, the gas drag force is given by the 
Epstein drag law 



-irp g s 2 VT^v 



(5) 



where p g is the gas density and vt is 4/3 times 
the mean thermal velocity. For the gas with mean 



2 Optical data arc available from http://www.astro.uni- 
jena.de/Group/Subgroups/Labor/Labor/labor.html 



molecular weight p g and temperature T, vt is 

4 / 8kT x 1/2 
vt = - 



3 \7r/i g mH 



(6) 



where k is Boltzmann's constant and hih is the 
mass of the hydrogen atom. In the case of Av — 
\Av\ ^> vt, the gas drag force is 



F g = -np g s 2 AvAv . 



(7) 



A standard way to connect these two regimes is 
(Kwok 1975) 

F g - -TT Pg s 2 (vT + Av 2 ) 1/2 Av . (8) 

The stopping time by the gas drag force is 
defined as t s = mAv/\F g \. It can be made 
non-dimensional by writing t s = Tgil^ 1 , where 
fix = (GM/r 3 ) 1 ' 2 is the Keplerian angular veloc- 
ity. The non-dimensional stopping time can then 
be expressed as 



(9) 



VI + Av 2 /v T 



where the square root factor represents the super- 
sonic correction. The subsonic stopping time pa- 
rameter equals 



4pdSVK 

3p g rvT 



(10) 



where vk = (GMjr) 1 ! 2 is the Keplerian veloc- 
ity. When T ss = 1, particles move slowly through 
the gas and couple to it in one dynamical time 
scale. Such particles are said to be marginally 
coupled to gas via gas drag, as opposed to the 
T ss 3> 1 case, representing large grains, decou- 
pled from the gas disk. Notice that to within 
a factor close to unity p g rvT/vK equals the gas 
surface density S g in a disk. This allows us to 
express T ss as T ss s=a Ei p /E g , the inverse ratio 
of Eg and a "surface density of one solid parti- 
cle", Sip = 4pds/3 (mass/cross-sectional area for 
a spherical particle). For example, a hypotheti- 
cal uniform surface density gas disk with radius 
100 AU and total mass 3 x 10" 5 M Q , or 10 Earth 
masses, would have S g = 8.5 x 10~ 3 g/cm 2 , the 
same as the Si p of an s = 51 p,m particle with den- 
sity pd = 1-25 g cm~ 3 . All the dust grains smaller 
than about 51 pm would then be well coupled to 
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gas in such a disk, while all sand grains would have 
stopping times longer than one orbital period. 

The velocity of gas v g can be written using the 
Keplerian velocity vk as 



v g =v K {l-rj) 



1/2 



(11) 



Quantity rj is the ratio of the pressure gradient 
force to the gravitational force, equal to 



n = - 



ril^pg dr 



(12) 



where r is the distance from the central star, and 
P g is the gas pressure. If gas pressure decreases 
with radius, which normally happens everywhere 
in the disk except at the inner disk edge, r\ > and 
the gas orbital motion is slightly sub-Keplerian, 
that is slower than Vk- 

In this paper, we consider the dust particles 
near the mid-plane of the gas disk, where we in- 
troduce cylindrical coordinates (r, 9) , and neglect 
the fact that the grains might have finite orbital 
inclinations. This is equivalent to considering the 
vertical motions to be averaged out by vertical in- 
tegration over an ensemble of particles in a short 
column above and below the point (r, 9). This ap- 
proach is justified by the generally weak coupling 
of vertical and planar motions of a grain in a thin 
disk. The equations of motion of a dust grain at 
r = (r, 0) with velocity v — (v r ,vg) are 



~dt 



r 



1-/3 1 



2v, 
c 



V K V r 



d , \ $v\v e v K 

dt {rV9) = —-i;^-^ 



, (13) 



(14) 



where T s given by equation (9) depends on the 
velocity of the grain. 

3. Radial Migration of Solids 

Small solid bodies migrate in disks with a speed 
depending on their size, radiative support and gas 
drag. If the radiation pressure exceeds central 
star's gravity (/3 > 1), then the grains (called (3- 
meteoroids) are blown away on a dynamical time 
scale, unless strongly braked by friction against 
gas. We do not concern ourselves with the observ- 
ability of rapidly escaping /3-meteoroids in this pa- 
per. (We do take into account their collisions with 



stable disk particles, however.) In this section, we 
consider grains that have nearly circular orbits and 
spiral in (or out) because of the drag forces. 

If the angular velocities of grains are larger than 
that of gas, the grains experience headwind and 
lose angular momentum to move inward. Inward 
migration of a grain subject to gas drag, but no 
radiation pressure or the PR drag, was discussed 
by Adachi, Hayashi, & Nakazawa (1976) and Wei- 
denschilling (1977). Exposure of dust grains to the 
radiation pressure often causes the orbital speed 
of grains to be slower than that of gas. In this 
case which we study below, the back- wind gives 
the grains angular momentum so that grains move 
outward. Below, we study several regimes of mi- 
gration allowing analytical solutions of the equa- 
tions of motion. 

3.1. Strongly Coupled Case (T s < 1) 

Consider the limiting case T s <C 1, in which the 
stopping time t s is much smaller than the dynam- 
ical time scale 0^ . As will be clear from the re- 
sults, the particle moves azimuthally at nearly the 
gas velocity, so that vq m v g . The grain assumes 
a finite terminal radial velocity v r when the radial 
forces come into balance so that dv r /dt = 0. The 
terminal velocity follows from equation (13), 



2 2 



VKVr 



0. (15) 



The term 2v r /c is negligible. From equations (9) 
and (11), the above equation becomes 



1 / v^ 1/2 

Solving this equation, we obtain 
v r = sgn(/3-ry)-^= 



2T SS ((3-V) — 







1/2 



(16) 

(17) 

1/2 



The particle migrates inward if /3 < r], or outward 
if > rj. 

The velocity of the grain in the azimuthal direc- 
tion vg is almost equal to the gas velocity v g , but 
there is a slight difference between these velocity, 
Ave = vq — v g . The angular momentum of the 
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grain is approximated by H = mrv g . The radial 
motion changes the angular momentum secularly 
at a rate 

dH mv g v r ( r drf 
dt 



(18) 



2 V 1 - V dr , 
The torque due to the gas and PR drag forces 
equals 



dH 
~dt 



m 

Tss 



1 + 



Avi 



1/2 



v K Ave-m0 c v K v s 



(19) 

where C — 0vk/c. As seen from equation (21) 
below, Ave <C v r . Thus, Avg in the parenthesis 
in the above equation can be neglected. Equating 
the above two expressions for angular momentum 
change and solving for Ave, the azimuthal wind 
velocity becomes 



Ave = -T aa [l 



-1/2 



(20) 



v r 

2v~^ 



1 - 



dr\ 



1 — r\ dr 



+ Pc 



where v r is given by equation (18). Notice that 
the formulae in this subsection remain valid in the 
supersonic migration regime, such as might arise 
from tight coupling of dust grain to gas, but even 
stronger coupling with stellar radiation field. 

3.2. Weakly Coupled Case (T s > 1) 

l{ > 1, grains are blown away quickly by 
strong radiation pressure force. In the case of 
S> 1, the radial velocity of a grain was derived 
by Lecavelier des Etangs, Vidal-Madjar, & Ferlet 
(1998). In this subsection, we consider the motion 
of a grain with < 1. The grain rotates with the 
speed ve = i>k(1 — P) 1 ^ 2 and angular momentum 
H = mrvyi{l — 0) 1 / 2 . The velocity difference from 
the gas is Ave = v K [(l - p) 1 / 2 -'(1 - r?) 1 / 2 ]. The 
torque by the gas drag and PR drag forces is 



dH 
~dt 



Av 2 



1/2 



v K Ave-m0 c v K ve 



(21) 

Because v r <C Ave (sec cq. [22] below), v 2 in the 
right hand side is neglected. The torque causes 
radial migration with velocity 



1/2 



1 



I] 



1-P 



1/2 



1 



v K -20 c v K 
(22) 



If PR drag force is negligible (the second term in 
eq. [22] is neglected), a grain of < r\ migrates 
inward, and a grain of > r\ migrates outward. 

3.3. Connecting Formulae For Arbitrary 
Stopping Time 

If the motions of grains and gas are close to 
the Kcplerian circular motion, we can derive ap- 
proximate formulae connecting the two regimes of 
weak and strong coupling of grains with gas. So- 
lutions of this kind, clarifying the size-dependent 
nature of migration, were presented for the clas- 
sical inward migration of grains in optically thick 
disks by Adachi et al. (1976) and Weidenschilling 
(1977). 

We assume that 7/, Ave/vx, v r /vK, and C = 
Pvk/c are much smaller than unity. Based on 
the assumption that grains undergo a steady-state 
migration in a disk without executing large radial 
motions in one turn of the orbit, we replace the full 
time derivatives d/dt in the equations of motion 
(13) and (14) with the advective terms v r d/dr. In 
other words, we set d/dt = 0. Then, neglecting 
the second order terms of ij, Ave/vK, v v /vk, and 
P c , the azimuthal equation (14) becomes 



Ave = —T a 



v r 

T 



PcVk) 



(23) 



The radial equation (13) becomes 



v K {0-v) + 2 Av e -—v r = 0. (24) 

1 s 

Substituting equation (23) into equation (24), we 
obtain 

P- V - 2p c T s 
Vr = T s + IT 1 VK ' 
Note that the stopping time T s defined by equa- 
tion (9) is a function of v r and Ave- Thus, v r 
and Ave cannot be expressed analytically. If the 
velocity difference of a grain with respect to gas 
is subsonic, T s « T ss is independent of the ve- 
locity of the grain. In this case, v r and Ave can 
be calculated analytically by replacing T s by T ss in 
equations (23) and (25). In the supersonic regime, 
v r and Ave must be found numerically by an iter- 
ative correction method. 

3.4. Steady Orbits of Grains 

At some distance from the central star, the 
grain may have a steady circular orbit without 
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radial migration or any torque. On such an or- 
bit, dv r /dt = v r = and dvg/dt = 0. From the 
r-component of equations of motion (13), the az- 
imuthal velocity of the grain is 

v e = v K (l - /3) 1 / 2 . (26) 

From the ^-component of equations of motion 
(14), 

(v e - v e ) + l3 c T s v g = . (27) 

The radius of the steady orbit is calculated as a 
solution of equation (27) for a given size of a grain 
s (or (3). If PR drag force is negligible ((3 C T S <C 1), 
equation (27) yields, vg — v s , and consequently 

= V- (28) 

Thus, a grain has the steady circular orbit at such 
a distance r from the star, where (3 — f](r). 

3.5. Size Dependence of Dust Particle Ve- 
locity 

In a gas disk, small grains of (3 > r] circu- 
late faster than the gas and migrate outward, 
while large grains of (3 < r\ rotate slower than 
gas and migrate inward, if PR drag force is neg- 
ligible. Figure 2 shows the dependence of the 
velocity of grains on their size. In this figure, 
p g = 2.18 x 1CT 17 g cm" 3 , v T = 1.07 x 10 5 cm s" 1 , 
and 77 = 2.24xl0~ 2 . These are the values the disk 
in model 1, described below in §4.1, has at 70 AU. 
The radial velocities v r are shown by circles. The 
filled circles represent that grains move outward 
(v r > 0), while open circles represent that grains 
move inward (v r < 0). These symbols are plot- 
ted by monitoring the velocities of grains in the 
numerical calculations of their orbits, described in 
§5. The grains smaller than 174/im migrate out- 
ward, while larger grains migrate inward. The ve- 
locity differences Avg in the azimuthal direction 
between grains and gas are shown by diamonds. 
The grains smaller than 174/xm rotate faster than 
gas, while larger grains rotate slower. The 174/im 
grains have (3 = 7], and thus steady circular orbits 
without radial migration. 

The grains with the size 0.1 < s < 10/im have 
similar values of the radial velocity. These grains' 
(3 are much larger than rj and inversely propor- 
tional to the grain size s (see Fig. 1). The stop- 
ping time of these grains are much smaller than 



the orbital period (T s — 1 at s « 30/im) and pro- 
portional to s (eq. [10]). It is seen that the radial 
velocity is constant in this size range, following 
from equation (25). 

The solid lines show the approximate formulae 
(23) and (25). The plotted values were obtained 
numerically after several iterations for the super- 
sonic correction. The dashed lines show the first 
guesses of iteration, which are calculated using the 
subsonic stopping time T ss instead of T s . The de- 
viations of the approximate and exact formulae 
for the 0.1 < s < lOyum grains are due to the fact 
that their radial velocities are comparable to the 
thermal speed of gas. 

4. Segregation of Dust According to Size 

As shown in §3, dust grains migrate inward or 
outward until they feel no net torque. The radius 
of the steady orbit is the function of grain's size 
s, which is obtained as the solution of equation 
(27). Thus, in a gas disk, grains segregate accord- 
ing to their sizes. In this section, we investigate 
the segregation of grains in a disk with fixed gas 
profile 3 . 

4.1. Models of Gas Disks 

We assume that the temperature and density 
structures of the gas disk have power law forms. 
The temperature structure of the disk that is 
transparent to the central star's visible radiation 
is (Hayashi 1981) 




where power law index q = —1/2. The mid-planc 
density has also power law, and at the outer radius 
of the disk r ou t , the disk has an edge of characteris- 
tic width Ar out . The density structure is modeled 
by 

» = »(T7u)'K t " h: sfr + 1 )- (30) 

3 Gas profile will change significantly in time if the total gas 
mass is much smaller than the dust mass, a case which we 
do not consider here. Gas migration and dust migration 
have opposite directions and magnitudes following from the 
conservation of total angular momentum. 
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In this model, the sound and thermal velocity c s 
and vt are, respectively, 



/ TkT \ 1/2 



Cs = 



V Mg TO H 



) 



= 1.17 x 10 5 



T 2.34 / L 



1.4 Hg \L & 



1/4 



1/2 



(r 



— j cms- 



(31) 



«T = 



8kT 
3 V 7r A'g TO H 

2.11 x 10 5 



1/2 



2.34 / L 



L 



1/4 



-,1/2 



1 AU 



9/2 



cm s 



(32) 



where T is the adiabatic exponent. The ratio of the 
pressure gradient force to the gravitational force rj 
is (cq. [12]) 



V = 



p + q- 



r sech x 



Ar out (tanhx + 1) 



(33) 

where x = (r out - r)/Ar out . 

In model 1, we adopt following assumptions. 
The disk gas consists of molecular hydrogen and 
atomic helium (p g = 2.34 and T = 1.4). The 
central star has the mass and luminosity of HR 
4796A: M = 2.5M Q and L = 21.0X©, respec- 
tively. The power law indexes of the density 
and temperature are p = —2.25 [which corre- 
sponds to the power law index of surface density, 
p s = P + (q + 3)/2 = -1.0] and q = -1/2, re- 
spectively. The outer radius of the gas disk is 
Pout = 100 AU. The width of the edge of gas 
disk is Ar out = C ou t(c s A^K)r 0Ut r ut- The crite- 
rion for the Rayleigh stability at the edge of the 
disk requires C ou t ^ 1- We adopt C out = 1.05, 
which yields Ar ou t = 12.1 AU. The mass of 
the gas disk is M g = 10M® that corresponds to 
po = 3.12 x 10~ 13 g cm~ 3 . The dust grains have 
material density of pd = 1.25 g cm -3 . We also 
consider the other models of the gas disk: smaller 
gas mass (M g = lM ffi ; model 2), steeper and gen- 
tler density gradient (p = —2.75, p s = —1.5; model 
3, and p = —1.25, p s = 0; model 4), and wider 
edge of the disk (C ou t = 2.0; model 5). The values 
of model parameters are summarized in Table 1. 



4.2. Segregation of Grains 

Figure 3 shows the ratio of pressure gradient 
force to the gravitational force r\ as a function of 
radius r. The solid line shows r\ for model 1. For 
comparison, rj of the gas disk with a wider edge 
(model 5) is shown as the dashed line. The value of 
the gas mass does not affect r\ and model 2 (M g = 
1M 9 ) has the same profile of r\ as that in model 1. 
The profiles of rj in model 3 and 4 (different density 
gradient; p = —2.75 and —1.25) are similar to that 
in model 1. In model 1 (solid line), r/ increases 
gradually with radius to 80 AU, and then more 
rapidly at the edge of the disk around 100 AU, 
because of the steep decrease in the gas density 
and pressure there. A wider disk edge (cf, the 
dashed line) causes a more gradual rise of rj. 

The radius of the steady orbits of grains are 
calculated from equation (27) and are plotted in 
Figure 4 (we plot the grain size s as a function 
of radius r). In our models, the PR drag force is 
negligible compared to the gas drag force except 
when the grains are at the outside of the disk, 
r > 130 AU. The dotted line shows the case in 
which the PR drag force is not included in model 
1 and coincides with solid line for r < 130 AU. 
Thus, the PR drag force can be neglected and the 
grains segregate to satisfy — rj. The smaller 
grains stay at farther from the central star. At the 
edge of the disk where rj varies rapidly (80 < r < 
110 AU), grains in the size range 10 < s < lOO^m 
find stable orbits. Again, a wider edge in model 5 
(dashed line) gives rise to a milder concentration 
of grains. 

5. Orbital Evolution of Dust Grains 

In §3, we derived the radial velocities of grains 
analytically assuming the orbits of the grains are 
nearly circular. In this section, we calculate the 
orbital evolution of the grains numerically by in- 
tegrating the equations of motion (13) and (14). 
It will be shown that the grains migrate to settle 
down in their own steady circular orbits, except 
at the outer edge of the gas disk, where the eccen- 
tricities of the grains' orbits are pumped up and 
dust disk extends beyond the gas disk. 
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5.1. Evolution of Orbital Radii and Eccen- 
tricities of Grains 

In Vega-type and transitional disks, dust grains 
are produced as ejecta of collisions between much 
larger parent bodies such as planetesimals, comets, 
or meteorites. These parent bodies are not af- 
fected by radiation pressure and have Keplerian 
orbits. If the velocity of ejection is much smaller 
than the orbital speed, the grains produced by col- 
lisions have the same Keplerian velocities as their 
parent bodies have. Because grains experience ra- 
diation pressure, the velocities of the grains are 
larger than that the grain on the circular orbit 
would have. Thus, the grains are endowed with 
orbital eccentricities e w (3/(1 — (3) (the equal- 
ity becoming exact, if parent orbits are circular). 
However, if the density of gas is high enough that 
the stopping time is much smaller than the or- 
bital period, the eccentricities of grains decrease 
quickly by the gas drag force and grains' orbits 
are circularized. In the inner part of our model 
disks (r < 20 AU), the orbits of grains with size 
s < 100/im are circularized rapidly (T ss < 1). 
Then, the grains spiral in or out as discussed in 
§3- 

Figure 5 shows the evolution of the distance of 
grains from the central star, r. The grains are 
assumed to be produced at 10 AU with eccentric- 
ity e = /3/(l — (3). These initial values are not 
important because the grains rapidly forget their 
initial states due to the gas drag force. The or- 
bits of the grains are calculated by integrating 
the equations of motion (13) and (14) numeri- 
cally by 7th-order Runge-Kutta method. The in- 
tegration are performed until t = 1000P ou t, where 
-Pout = 27r/QK(?*out) = 632yr is the orbital pe- 
riod at r out (= 100 AU). Figure 5a shows the 
evolution of orbits in model 1 . The orbits of the 
grains are quickly circularized and the grains mi- 
grate outward. The grains larger than 100/mi ar- 
rive at their steady orbits in ~ 100P ou t, and settle 
down in these orbits. Smaller grains migrate to 
farther locations. A grain of 20^m arrives at the 
radius of the steady orbit that is in the outer edge 
of the disk in lP ou t, then its orbital eccentricity 
is pumped up. The excitation of the eccentric- 
ity is seen in Figure 5a as the oscillation of the 
distance r. Finally, the orbital eccentricity of the 
20yum grain decreases due to the gas drag force, 
and the grain stays in the circular steady orbit. 



The orbital eccentricity of a grain of 8/im is also 
pumped up when the grain arrives at the edge of 
the disk, and then decreases. Because the 8/itm 
grain spends most of its time staying outside the 
disk (it is in the gas disk only when being just 
around the pericenter of the orbit), however, the 
gas drag force is not strong enough to damp the 
eccentricity quickly. Thus, the 8/im grain stays in 
eccentric orbit for long time (more than 1000P ou t)- 
The evolution of orbits in the disk with smaller gas 
mass (M g = 1M ; model 2) is shown in Figure 
56. The main feature of the evolution is similar to 
that in model 1. The settling time for grains to 
arrive at the steady orbits is longer than that in 
model 1; ~ 100P ou t for 100 and 300/xm grains and 
- 1000P ou t for 200/im grains. 

5.2. Excitation of Eccentricity at the 
Disk's Outer Edge 

The orbital eccentricities of grains smaller than 
20/im are excited when these grains arrive at the 
outer edge of the gas disk. These grains have large 
(3 (0.19 for 20/zm and 0.47 for 8/im grains). Un- 
til the grains arrive at the disk edge, the gas drag 
force on such the small grains are strong enough 
that the grains are trapped by the gas and move 
with the gas. The gas rotates faster than the 
grains (t] = 0.10 at r out ), and the grains gain 
angular momentum from the gas and move out- 
ward. Because the gas density decreases rapidly 
at the disk edge, the gas drag force ceases sud- 
denly when the grains pass through the disk edge. 
Therefore, the grains are launched into the Kepler 
orbit (modified by the radiation pressure) with ini- 
tial velocity of v g = (0,Wk[1 — v] 1 ^ 2 ) & t the disk 
edge. With this initial velocity, the orbits of the 
grains are eccentric with e = ((3 — f])/(l — [3) and 
pericenters at the disk edge. After launching into 
the eccentric orbits, the grains experience the gas 
drag force that works to circularize the grains' or- 
bits only when they are around their pericenter 
at the disk edge. Thus, the eccentricities of the 
grains decreases slowly. 

In the above discussion, we assumed that the 
gas density in the disk is high enough for the gas 
drag force to hold the grains tightly with the gas. 
We arrive at the same conclusion when the gas 
density is low and the gas drag force is not strong 
enough to hold the grains. Consider a small grain 
whose (3 is much larger than r\ at the disk edge, 
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and assume that the grain has slightly eccentric 
orbit. The semi-major axis and eccentricity of its 
orbit are a and e, respectively. Because (5 3> f], 
the grain rotates much slower than the gas, and 
experiences back- wind throughout its orbit, even 
when it is at the pericenter. The evolution of the 
eccentricity of the orbit by the perturbating drag 
force is given by the Gauss's equation (Brouwer & 
Clemence 1961) 



de 



1 



1-e 1 



1/2 



uk (a) V 1 - P 
x F r sin ip + F^ 



(34) 



cos ip + e 
1 + e cos ip 



+ cos ip 



where ip is true anomaly, and F r and F^ are the 
radial and azimuthal components of gas drag force, 
respectively. At the pericenter (ip = 0) and the 
apocenter (ip = it), equation (35) becomes 



de _ 2 
dt vk(o) 



l-P 



1/2 



(35) 



where plus and minus signs correspond to the peri- 
center and apocenter, respectively. The grain ex- 
periences back-wind throughout the orbit. Thus, 
F^/j is positive both at the pericenter and apocen- 
ter. From the equation (35), it is seen that the ec- 
centricity increases at the pericenter and decreases 
at the apocenter. When the grain is in the edge of 
the disk, the gas density at the pericenter is much 
higher than that at the apocenter. Therefore, the 
grain experiences stronger back-wind at the peri- 
center and its eccentricity increases. The increase 
in the eccentricity accompanies the increase in the 
velocity at the pericenter and continues until the 
velocity arrives at the same velocity of the gas. 
When the velocity of the grain at the pericenter 
becomes the same as that of the gas, the eccen- 
tricity has grown up to bee = (P — rf) / (1 — P) . Af- 
ter the eccentricity arrives at the maximum value, 
it starts to decrease. This is because the gas 
drag force at the pericenter pushes the position 
of the pericenter slightly outward. The eccentric- 
ity evolves to keep the velocity of the grain at the 
pericenter being equal to the velocity of the gas 
and maintains the value of e = (f}—rj)/(l — (i). Be- 
cause r] is an increasing function of r, the value of r\ 
at the pericenter increases as the pericenter moves 
outward, and the eccentricity decreases. The gas 
density at the pericenter and the gas drag force 



decrease as the pericenter moves outward, and the 
speed of the eccentricity damping becomes more 
slowly. 

Figure 6a illustrates the evolution of the eccen- 
tricity and the semi-major axis of an 8/zm grain. 
The average of the change in the eccentricity over 
one orbital period P is calculated by 



-- f 



de 
~dt 



dt 



(36) 



The integration is performed numerically and the 
effect of the PR drag force is included. The 
Gauss's equation of the variation in the semi- 
major axis is 



da 
~dl 



fi K (a)(l-/3) 1/2 (l-e 2 ) 1 /2 



F r es\nip + F$ 



a(l-e 2 ) 



■ (37) 



The average over one orbital period, (da/dt), 
is calculated in analogy with equation (36). The 
evolutions of a and e are shown in Figure 6a as 
arrows. The length of the arrows are set to be 
constant and does not represent the rate of the 
variation. A grain inside the disk (a < 100AU) 
moves outward and its eccentricity decreases. If 
the grain's pericenter is on the edge of the gas 
disk (100 < a < 200AU and < e < 0.6), its ec- 
centricity grows. For a grain whose pericenter is 
outside the disk, the eccentricity and semi-major 
axis decrease. The evolution of the orbit of a grain 
whose initial a and e are 100 AU and 0.05, respec- 
tively, is plotted as circles on Figure 6a. The time 
interval of plotting circles is logarithmic: the first 
time interval is 0.01P O ut and each time interval is 
1.2 times the previous one. The last circle repre- 
sents the orbit at t = 1000P O ut- First, the grain's 
semi-major axis and eccentricity grow, then the 
eccentricity turns to decrease. 

In Figure 66, the evolution of orbital elements 
of a 100/xm grain is shown for comparison. It is 
seen that the eccentricity always decreases and the 
semi-major axis changes toward the value of the 
steady orbit. Circles show the evolutional path of 
a grain with the initial semi-major axis a = 100AU 
and eccentricity e = 0.5. The eccentricity de- 
creases quickly, and the semi-major axis decrease 
at first, and then increases as the grain approaches 
the steady orbit. 
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5.3. Radial Distribution of Dust Grains 

The radial distribution of dust grains are ob- 
tained by the numerical calculation described in 
§5.1. Figure 7 shows the positions of the peri- 
center and apocenter at t = 1000P ou t against the 
grain size. The solid line corresponds to model 
1. The grains smaller than Q.7/im are blown away 
(launched into hyperbolic orbits when they move 
to outer edge of the gas disk). The grains of 
6.7 < s < 12^m have eccentric orbits (e > 0.01) 
with their pcriccntcrs on the outer edge of the gas 
disk. The grains larger than 12^im have circular 
orbits. The radii of the circular orbits coincide to 
the radii calculated from the condition that the 
orbits are steady (eq. [27]) and shown as the dot- 
ted line. The dashed line corresponds to model 2 
(smaller gas mass; M g = lM ffi ). Because the gas 
drag force is smaller than that in model 1, the ec- 
centricities of small grains remain larger. The size 
of the smallest grain having circular orbit is 13fim, 
slightly larger than that in model 1. 

In summary, the migration of dust grains pro- 
duced at the inner part of the gas disk proceeds 
as follows. First, the grains migrate outward. The 
large grains of s > 12^m find the steady circu- 
lar orbits in the gas disk and settle down. When 
the small grains (s < 12^m) arrive at the edge of 
the gas disk, their eccentricities are excited. The 
grains of 6.7 < s < 12/im have eccentric orbits, 
and make dust disk extended outside the gas disk. 
The grains small enough for their eccentricities to 
be pumped up to e > 1 are blown away. The grains 
of 10 < s < lOO^m are concentrated at the edge 
of the gas disk. It is expected that these grains 
make a dust ring at the edge of the gas disk. In 
the next section, we estimate the number density 
of grains in the ring. 

6. Lifetime and Density of Dust Grains 

The grains larger than about 10/im have steady 
circular orbits in the gas disk and segregate ac- 
cording to their sizes. The size of grains in the 
steady orbit is function of distance from the star 
r, written as 

s s = s s (r) . (38) 

The grains of size s s in the steady orbit do not 
collide each other, because their orbits are circu- 
lar. However, smaller grains (s < s s ) whose steady 
orbits arc farther from the star flow outward and 



cross the orbit of grains of size s s . These smaller 
grains collide and destroy the grains in the steady 
orbit. In this section, we calculate the lifetime of 
grains in the steady circular orbit in the gas disk, 
then estimate their number density. The small 
grains (s < 10^m) extended outside the gas disk 
have eccentric orbits and collide each other. The 
calculation of the lifetime of these small grains re- 
mains as future work. 

The grains arc originally produced by collisions 
of parent bodies such as planetesimals or comets. 
In this section, we assume that these parent bod- 
ies are distributed in the innermost part of the gas 
disk (say, <~ 10AU). The grains produced by colli- 
sions migrate outward to settle down in the steady 
orbits. When the grain on the steady orbit is bro- 
ken up by collisions, the fragments flow outward 
and settle down in the new steady orbits. Thus, 
grains gradually migrate outward every collisional 
destructions. The steady orbit of grains of size 
10 < s s < 100/zm are concentrated at the edge of 
the gas disk (Fig. 4). Thus, it is expected that 
the grains of 10 < s s < 100^m gather in the edge 
of the gas disk and make a dust ring. The grains 
are finally blown away when collisional destruction 
reduces their sizes so small that they cannot have 
steady orbit in the gas disk. 

6.1. Collisional Destruction of Grains 

There are two phases of collisional destruction. 
A grain is broken up by a collision with large 
impact kinetic energy (catastrophic disruption), 
while part of a grain is ejected by a collision with 
small energy (cratering) . Consider the destruction 
of a grain of mass m — 4irpdS 3 /3 by the collision 
with a grain of mass m! . The mass Am lost from 
the target grain m by cratering with the collisional 
velocity u is (Holsapple 1993) 




where K is a coefficient representing efficiency of 
cratering and the value of the exponent fj, is un- 
certain now. We assume that the mass loss Am 
is proportional to the impact energy and adopt 
/i = 2/3. In equation (39), Am/m becomes larger 
than 1 for large u, that corresponds to catastrophic 
disruption. Of course, Am/m cannot be larger 
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than 1 and should be 

= max K— , I) . 40 

m I m \ cm s 1 / J 

The value of K is calculated from the minimum 
energy for the catastrophic disruption. Let Q* be 
the threshold of the specific energy per unit target 
mass for the catastrophic disruption. The thresh- 
old energy Q* depends on the impact velocity u 
and target size s, and written as (Housen & Hol- 
sapple 1999) 

/ a \ 2— 3ju / g \ <5 / g \ 5 

°* =0 Hkn7^J llOW = Q HlO^J 

(41) 

where we used [i = 2/3. The values of Qq are 
obtained by the experiments: Qq = 10 7 erg g _1 
for rocks (Housen & Holsapple 1999) and Q* Q = 

4 x 10 5 erg g _1 for ice (Fig. 9 in Arakawa 1999). 
The value of the exponent S is uncertain. Housen 
k Holsapple (1990) proposed that 6 = -0.24 and 
Housen & Holsapple (1999) recommended that 

5 = —0.67. We are interested in the lifetime and 
density of grains constituting the dust ring at the 
edge of the gas disk. The sizes of these grains are 
10 < s < lOO^m. Thus, we estimate Q* for the 
grain of size s = 50^m as 




6.2 x 10 7 — 1.6 x 10 9 erg g 1 for rocks 
2.5 x 10 6 - 6.5 x 10 7 erg g" 1 for ice 



(42) 

where smaller values correspond to 6 — —0.24 and 
larger values correspond to S = —0.67. The defi- 
nition of the catastrophic disruption is Am/m > 
1/2, which yields 

R = _J_ f 1.5 x 10~ 10 - 4.0 x 10~ 9 for rocks 
4Q* \ 3.8 x 10~ 9 - 1.0 x 10~ 7 for ice 

(43) 

6.2. Lifetime of Grains 

In this subsection, we derive the equation for 
the life time of a grain of size s s on the steady or- 
bit. A grain s s collide with smaller grains crossing 
its orbit. The probability of collision with grains 
in the size range [s', s' + ds'] in 1 sec is 

Pcoi = crun'ds' , (44) 

where a = ir(s s + s') 2 is the collisional cross sec- 
tion, n'ds' is the number density of grains in the 



size range [s',s' + ds'}, and u = \v ' — v\ is the 
relative velocity of the grains. The velocity of the 
grain s s is the same as the velocity of the gas, 
i.e., v = (0,t> g ). The velocity v' — (v' r ,v' ) of 
the grain s' is obtained by the orbital calculations 
of grains described in §5.1. For small grains in 
the inner part of the disk where the gas density 
is high, the non-dimensional stopping time T s is 
much smaller than unity. We use the analytical 
formulae (18) and (21) to calculate the velocities 
of small grains if T s <C 1, because numerical cal- 
culations of orbits consume much time. The den- 
sity n'ds' is calculated from the production rate 
of dust grains. Let the number of grains that are 
in the size range \s', s' + ds'} and cross the circle 
of the radius r (the orbit of the grain s s ) in 1 sec 
be Ap ro (r, s')ds' . The number density N plo (r,s') 
is assumed to have power law form in s'. We also 
assume that the grains are originally produced in- 
nermost part of the gas disk and flow outward. 
Only the grains smaller than s s cross the orbit of 
the grain s s . The grains larger than s s stay in 
their steady orbits which are closer to the star, 
but when they are broken to pieces smaller than 
s s , these pieces flow outward and cross the orbit 
of the grain s s . Thus, the grains crossing the orbit 
of the grain s s are constituted of smaller grains, 
and A^ pro (r, s')ds' are written as 

AT pro (r )S ')^'=(^ o(r) (-) a<is/ fors ^ Ss . 

( for s > s s 

(45) 

The mass of grains crossing the orbit of the grain 
s s in 1 sec is 

f Ss 4tt • / s' \ a 
M pro = J —s' 3 Pd N (_j ds' , (46) 

where s m ; n is the minimum size of grains. In the 
steady state, M pro is independent on r and repre- 
sents the mass production rate of grains in the gas 
disk. From equation (46), 

Note that N is the function of r because s s is the 
function of r. The other expression of number flux 
is iV pro (r, s')ds' = 2irrhdv' r n' ds' , where hd is the 
thickness of the dust disk. Solving this equation 
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and equation (45) gives 
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The mass loss rate from the grain s s with mass m 
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6.3. Spatial Density of Grains 



Consider an annular part of the disk between r 
and r + dr. The annulus is made of the grains of 
the size [s s , s s + (ds s /dr)dr]. The total number of 
the grains in the annulus is 



N d ds s — 2irrhdn s dr 



(51) 



where n s is the number density of the grains in the 
steady orbits 4 . The number of grains destroyed by 
collisions in 1 sec is 



N dcs ds s = 



N d ds s 2irrh d n s dr 
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t\if e 



(52) 



The supply of the grains to the annulus is 
N pro {r, s s )ds s = No{s s /cm.) a ds s - In the steady 
state, the rate of destruction and supply balance 
each other. From N dcs ds s = N pro (r, s s )ds s , we 
obtain 
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4 Note that the unit of n B is cm 3 , while that of n' is cm 4 . 



We define the optical depth of the dust disk in 
the vertical direction using the geometrical cross 
section as 



T = 



irs 2 h d n s 
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Note that the density and optical depth of the 
grains are independent on N or the production 
rate of grains M pro - 

6.4. Results of Numerical Calculations 

In this subsection, we calculate the life time of 
grains and density (optical depth) of the dust disk 
in model 1. We assume that the grains are made 
of rocky material and adopt K = 1 x 10~ 9 . The 
mass production rate of grains are uncertain. The 
models of dust disks (however, without gas disks) 
suggested that the rate of mass loss from the disk 
is of order 10 -6 M®yr -1 for (3 Pic (Artymowicz 
& Clampin 1997), and HR 4796 A (Augereau et 
al. 1999b). In the steady state, the mass loss 
rate balances with the mass production rate. We 
adopt mass production rate M pro = lO^Meyr -1 . 
The density of grains or the optical depth of the 
disk do not depend on the mass production rate, 
while the life time of grains is inversely propor- 
tional to it. The thickness of the dust disk is as- 
sumed to be h d = O.lr and the minimum size of 
grains s m j n = O.Ol^m is adopted. The power law 
index of A pro in the equation (45) is assumed to 
be a = —3.5, a typical value usually assumed for 
collisional cascade. 

The solid line in Figure 8a shows the life time 
of grains in model 1 (gas mass M g = 10M e ). The 
grains larger than 100/^m are in the inner part 
of the gas disk. Their life time is order of 1000 
times the orbital period P ou t at r out (— 100 AU) 
and decreases with the size because larger grains 
inhabit inner part of the disk where the density of 
colliding grains is higher (eq. [48]). The grains less 
than lOO^m are in the edge of the gas disk. At the 
edge of the disk, the collisional velocity of grains 
becomes large because the gas density is not high 
enough for the gas drag force to hold the grains 
in the same velocity. Thus, the life time of these 
small grains is much shorter than that of larger 
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grains (s > 100/im) in the gas disk. In model 
1, the life time of grains is much larger than the 
migration time to the steady orbits (the dashed 
line in Fig. 8a). Thus, the grains spend most of 
their time in their steady orbit and the dust disk 
is mainly composed of the grains in the steady 
orbit. The contribution of grains flowing toward 
the steady orbit to the density of the dust disk or 
the optical depth can be neglected. 

The solid line in Figure 9a shows the optical 
depth of the dust disk in model 1. The dominant 
profile of the optical depth is the peak at the edge 
of the gas disk (~ 90AU). This peak is made of 
the grains of 10 < s s < 100/im. Because these 
grains gather in the edge of the gas disk (see Fig. 
4), the density of grains becomes high and a dust 
ring forms there. The life time of grains outside 
the gas disk (r > 100AU and s < 10/im) is much 
shorter than that of grains in the gas disk. Thus, 
the optical depth declines steeply outside the gas 
disk. The large grains (s s < 300/^m) condense at 
the innermost part of the gas disk. The density of 
gas is high at the inner part of the disk, and strong 
gas drag force suppresses the collisional disruption 
of grains. Thus, the optical depth at small radius 
becomes large and an inner dust disk forms at r < 
40AU. 

6.5. Models of Dust Distribution 

The structure of dust disks strongly depends on 
the properties of gas disks. In this subsection, we 
show how the structure of the dust disk varies with 
the variation of the gas disk. 

6.5.1. Small Gas Mass 

The solid line in Figure 86 shows the life time 
of the grains in the gas disk with smaller mass 
(model 2; M g = lM e ). The life time of the 
grains is much shorter than in model 1, because 
the weaker gas drag force results in high collisional 
velocity of grains. In model 2, the life time is com- 
parable or smaller than the time for the grains 
to migrate to their steady orbits (the dashed line 
in Fig. 86). This means that the dust disk is 
composed of both the grains in the steady orbits 
and the grains flowing outward to the steady or- 
bits. Some grains flowing outward are destroyed 
before arriving at the steady orbit. In this case, 
the calculation of the density of grains described 



in §6.3 is not appropriate. The ring structure at 
the edge of the gas disk would be less prominent, 
because some of grains making the ring would be 
destroyed before arriving at their steady orbit. To 
obtain the structure of the dust disk, the calcula- 
tion taking into account the destruction of grains 
on the way to the steady orbit must be needed. If 
production rate of grains is 10 times smaller and 
M pTO = 10 _7 M©yr _1 , the life time of grains is 10 
times longer, and our calculation of density would 
be appropriate. The optical depth in the model 
of M g = 1M® is shown as the dashed line in Fig- 
ure 9a, assuming that the life time of grains are 
much larger than the migration time. The profile 
of the optical depth is similar to that in the model 
of M g = 10M ffi , while its magnitude is about 10 
times smaller. At the peak of the optical depth 
in the dust ring (r w 90AU), r = 1.1 x 10 _1 for 
M g = 10M© and r = 1.2 x 10~ 2 for M g = lM e . 

6.5.2. Different Gas Profiles 

In Figure 96, we show the optical depth of 
the dust disks for various density gradients of 
gas disks. The dashed line shows the disk with 
steeper density gradient (the power law index of 
the gas density p — —2.75 corresponding to that 
of the surface density p s = —1.5; model 3). The 
mass of the gas disk is the same as in model 1 
(M g = IOM0). Because of the steeper density gra- 
dient, the gas density in the inner part of the disk 
is higher than that in model 1, while the gas den- 
sity at the edge of the disk is lower. The lower gas 
density at the edge of the disk causes the shorter 
life time of grains and the smaller optical depth of 
the dust ring. The grains in the inner part of the 
disk are strongly held by dense gas. The collisional 
velocity of these grains are small and collisional de- 
struction of grains is not effective. The grains in 
the inner part of the disk survive longer than the 
grains in model 1, and the optical depth increases 
with decreasing radius r. The inner optically thick 
disk may obscure the star from the outer part of 
the disk. In this case, the grains in the outer part 
of the disk do not experience the radiation pres- 
sure force and do not have steady orbits, so that 
the dust ring would disappear. 

The dotted line in Figure 96 shows the disk with 
the constant surface density (p — —1.25 and p s = 
0; model 4). The gas density in the inner part 
of the disk is much lower than that in the model 
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1, while the gas density at the edge of the disk 
is higher. The optical depth of the dust ring is 
larger than in the model 1 because of the higher 
gas density at the edge of the disk. In the gas 
disk with the constant surface density, even the 
grains in the innermost part of the disk arc not 
strongly held by gas and are destroyed in short 
time by collisions with high velocity. The number 
density of these grains is not as large as that in 
model 1 or 3. The inner disk is optically thin. 
In model 4, because the optical depth decreases 
with decreasing radius r, the inner disk does not 
obscure the star. 

In the previous models, we assumed that the 
gas disks have sharp outer edges (marginally 
Rayleigh stable at the edge where C ou t = 1-05). 
The sharpness of the edge is not known well. Re- 
cent images of the dust disk of GG Tau show a 
very sharp edge (Ar out /r out <~ 0.04; Guilloteau, 
Dutrey, & Simon 1999). However, the edge of gas 
disk may be smooth. The optical depth of the dust 
disk in the gas disk with smooth edge (C ou t = 2.0; 
model 5) is shown in Figure 9c. The concentra- 
tion of grains at the disk edge is moderate, and 
the ring feature is weak. 

6.5.3. Grains of Smaller Collisional Strength 

The value of the efficiency of collisional destruc- 
tion K is uncertain (eq. [43]). If the grains are 
stronger (weaker) , their life time is longer (shorter) 
and their density is higher (lower). As an exam- 
ple, we show the optical depth of the dust disk 
with K = 2 x 10~ 8 , which is 20 times larger than 
K in model 1 and corresponds to the efficiency of 
collisional destruction of icy grains (model 6), in 
Figure 9d. The life time of the grains is about 
20 times shorter than that in model 1, and the 
optical depth is also 20 times smaller. The icy 
grains are destroyed more easily than the rocky 
grains and are blown away by the radiation pres- 
sure more quickly. Therefore, the dust disk are 
mainly composed of rocky grains with little icy 
mantle. This is consistent with the properties of 
the grains derived through the spectral fitting of 
thermal emission from the HR 4796A dust ring by 
Augereau et al. (1999b). 



7. Morphology of the Ringed Systems 

We consider first the global quantities (mass 
and luminosity) of the dust, as well as gas, in the 
models, and compare them with some observed 
quantities. Next, we discuss the morphology of 
the model disks in the thermal and scattered light 
images. 

7.1. Mass and Luminosity of Model Dust 
Rings 

The mass of the dust disk M<± is given by 

8 /"'"out 

M d = -ir 2 p d / rh d sln s dr . (55) 

J Tin 

We set r; n = 50 AU and r out = 150AU to calculate 
the mass of the ring part of the disk. The mass of 
the dust ring in model 1 (gas mass M g = 10M e ) is 
9.0 x 10 _1 M e and the dust mass in model 2 (M g = 
lM e ) is 1.0 x 10 _1 M e . The gas to dust ratio is 
about 10. Thus, the gas-dust coupling affects the 
motion of the gas only slightly, at best. (If the 
dust grains settle down in the mid-plane of the 
disk, the dust density may exceed the gas density. 
Even in this case, the dust grains do not affect the 
gas motion strongly, because of the corotation of 
gas and dust in steady orbits.) 

The above estimate of the mass assumes that 
the dust ring is composed of grains smaller than 
s <~ 200^m. Augereau et al. (1999b) suggested 
the presence of meter-sized bodies based on the 
spectral energy distribution of the HR 4796A dust 
disk. The gas drag and radiation pressure force 
have little effect on such large bodies, so unlike 
for the dust component, the distribution of these 
bodies is not strongly tied to the distribution of 
gas. The presence of large meteorites and plan- 
ctcsimals (which do not fall under the category 
of dust) necessarily increases the total mass of the 
disk of solid matter above our estimate of the dust 
disk mass. Additional calculations are needed to 
model the long-term evolution of the large bodies 
in disks. 

The total luminosity of thermal emission and 
scattered light from the dust ring is 

Ld= 2 L J rin — dr =2 L J nn r dr - {56) 

The above equation assumes the extinction cross 
section to be geometrical. The luminosity of the 
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dust ring in model 1 is 1.7 x 10 2 L and that in 
model 2 is 2.0 x 10~ 3 L. 

The shape of the dust ring depends on the 
sharpness of the edge of the gas disk (Fig. 9c) . If 
the edge is more gradual, the ring is wider and its 
peak density is lower. However, the mass and lu- 
minosity of the ring do not strongly depend on the 
sharpness of the disk edge. In model 5, where the 
width of the edge is about 2 times wider than that 
in model 1, the mass and luminosity of the ring are 
M d = 7.2 x lO^Mg, and L d = 1.4 x 10~ 2 L, re- 
spectively, which are only about 20% smaller than 
those in model 1. 

We can compare the above numbers with the 
masses of gas and dust in transitional disks. HR 
4796A dust ring has been observed by thermal 
emission (Jayawardhana 1988; Koerner et al. 
1998; Telesco et al. 2000) and by scattered light 
(Schneider et al. 1999). The emission from the 
dust ring has a peak around 70 AU (Schneider 
et al. 1999; Telesco et al. 2000). The luminosi- 
ties of thermal emission and scattered light are 
L th ^ 5 x 10~ 3 L (Jura 1991) and L sc ~ 2 x 10~ 3 L 
(Schneider et al. 1999), respectively. The total lu- 
minosity of the ring is Ld — 7 x 10~ 3 L. A similar 
dust distribution with a peak at 72 AU and lumi- 
nosity id - 8 x 10~ 3 L is found in our model 7, 
where the gas disk has a smooth outer edge at 80 
AU and mass M g = 4M® . The gas mass in model 
7, M g = 4M®, is consistent with the observational 
upper limit M g < 7M® (Zuckerman et al. 1995; 
Greaves et al. 2000). 

The luminosities of the thermal emission and 
scattered light of HD 141569A are L th ^ 8x 10~ 3 L 
(Sylvester et al. 1996) and L sc ~ (2 - 4) x 10~ 3 L 
(Augereau et al. 1999a; Weinberger et al. 1999), 
respectively. The total luminosity of the dust is 
Ld — 1 x 10~ 2 L. The model which reproduces 
the dust ring with a peak density at 325 AU and 
luminosity Ld — 1 x 10~ 2 L has the gas disk with 
an edge at 370 AU and mass M g = 50M® (model 
8). The inferred gas mass is consistent with the 
mass M g ~ 20— 460M® derived by the observation 
(Zuckerman et al. 1995). 

Although the gas masses derived by our models 
are consistent with the observations, we must note 
that any exact match could be accidental. Our 
models have large uncertainties. In particular, the 
efficiency of the collisional destruction K has an 
order of magnitude of uncertainty (cq. [43]). If K 



of the actual grains were larger/smaller than the 
value we adopted (K = l.Ox 10~ 9 ), then the grains 
would be more/less easily destroyed and the mass 
and luminosity of dust in a given gas disk would 
be lower/higher. 

7.2. Model Images of Scattered and Ther- 
mal Radiation 

In Figure 10, we present the theoretical im- 
age of model 4 (constant gas surface density), in- 
clined to the line of sight by 60° from pole-on, 
in the near-IR scattered light (gray-scale bright- 
ness map at 1.1 /im wavelength), and the mid- 
IR thermal radiation (over-plotted contours of the 
18.2 /im flux). Our modeling is somewhat sim- 
plistic, in that the scattering function is assumed 
to be isotropic (which does not alter the general 
morphology, just removes any up-down asymme- 
try of the theoretical image). However, the tem- 
perature profile as a function of radius needs to be 
computed as accurately as possible for a reliable 
prediction of thermal IR flux distribution, so we 
compute it within Mie theory, from the assumed 
properties of grains (specified earlier in §2). 

The dust morphology in the model is qualita- 
tively similar to those inferred from observations of 
the transitional system HD 141569A. The 18.2/xm 
flux is centrally concentrated (like that imaged by 
Fisher et al. 2000), while the scattered light shows 
two separate components: the inner disk and the 
outer ring, with the brightest parts being the two 
projection-enhanced ansae. If these two compo- 
nents are largely overlapping, they may create the 
impression of a shallow gap, resembling qualita- 
tively the image at A « 1.1/um by Weinberger et 
al. (1999). [However, see Augereau et al. (1999a) 
for the dust distribution from their observations 
at A ~ 1.6/im.] The model disk in Figure 10 has 
a constant gas surface density with a downturn 
beyond 100 AU. Model disks with surface density 
decreasing with radius (model 1 etc.) show much 
more prominent inner disks both in scattered light 
and in thermal emission. Comparison with the ob- 
servations suggests that gas is not concentrated in 
the inner part of the disk. The flat profile of gas 
surface density might arise after the viscous evo- 
lution. 

In order to achieve a quantitative agreement of 
our models with the observations of HD 141569A, 
we would need to: (i) consider non-isotropic scat- 
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tering function, (ii) allow for a more gradual outer 
gas truncation, and (iii) convolve the theoretical 
images with the appropriate point-spread func- 
tions. As already mentioned, we defer the compre- 
hensive modeling of individual objects to a subse- 
quent paper. 

Another system which might be qualitatively 
described by our dust migration models, HR 
4796A, shows an enhanced-brightness ansae of 
the dust ring in scattered light (Schneider et al. 
1999), similar to our model in Figure 10. The 
18.2/im images also indicate a ring-like morphol- 
ogy (Jayawardhana et al. 1998; Koerner et al. 
1998; Telesco et al. 2000). In other words, the 
dust ring with characteristic radius of 70 AU, pro- 
ducing the scattered light, also produces most of 
the 18.2/im thermal radiation. 

Our thermal mid-IR map shows a more cen- 
trally concentrated morphology, therefore our 
modeling procedure needs to be argued in order to 
fit the observations. Two possibilities exist. First, 
our model image represents the emission from 
dust in steady circular orbits, and omit the possi- 
ble contribution by other particles, those flowing 
out on spiral orbits due to (3 > 1, and the ones 
pumped up to high eccentricity at the edge of the 
disk (see §5.2). The combined area of such dust 
grains is small, but they could in principle be seen 
because of the small physical sizes, and the associ- 
ated higher temperatures (see Fig. 4 in Augereau 
et al. 1999b and Fig. 5 in Wyatt et al. 1999). 
We note that the thermal emission in mid-IR (the 
Planck function in Wien regime) and thus the 
visibility of grains in the outer ring is often very 
sensitive to the dust grain temperature; it is not 
uncommon to find that the predicted 18.2/im flux 
from the ring becomes 6 times under just ~20% 
increase in dust temperature (from 70K to 84K), 
which could make the model ring visible. Another 
way to assure a quantitative fit to HR 4796A data 
is a modification of the modeled gas profile, which 
depletes the gas inside the ring radius. Our model 
would then produce a more sharply defined ring 
and less dust inside the ring, because, at the re- 
gion where gas is depleted, collisional disruption of 
dust grains is violent. This depletion of gas (and 
dust) might in this case be due to one or more 
planets or a low-mass companion inside about a 
half of the ring radius (~ 40AU). 



8. Discussion and Conclusions 
8.1. Dynamics of Dust Migration 

We have studied in the present paper the com- 
bined action of gas drag and radiation pressure 
on dust grains in optically thin circumstellar gas 
disks. Our results on dust migration are as follows. 

1. Grains can migrate inward or outward, de- 
pending on whether its azimuthal motion is faster 
or slower than the rotation velocity of gas. Even 
grains smaller than the blow-out radius {(3 > 1) 
have spiral trajectories in gas disks. 

2. We have provided analytical formulae for the 
radial migration speed in §3, which agree with the 
detailed numerical integrations. The radial migra- 
tion speed depends on gas density, but could be 
large. For grain sizes s < 100/im in disks around 
A-type stars of radius comparable with 100 AU 
and masses of order 1 — 10 Earth masses it reaches 
a significant fraction of gas sound speed. 

3. At some radius in the gas disk, a particle may 
reach a stable equilibrium orbit where it corotates 
with the gas. Equilibrium is always found outside 
the radius where the gas pressure is maximum. 

4. The radius of the steady orbit depends on 
the grain size. That dependence is given by the 
rotation curve of the gas disk, depending on the 
radial gas distribution. Smaller grains normally 
find a steady orbit farther from the central star. 
This results in the fractionation of grains accord- 
ing to their sizes. 

5. The eccentricity of a grain is strongly 
damped through out most of the disk, but is 
excited and maintained at a value up to e = 
([3 — — (3) at the edge of the gas disk, where 
the gas density falls rapidly (e.g., as a steep power 
law or exponentially). The rule number 3 above 
applies then, not to the whole orbit, but to the 
vicinity of the pericenter of the orbit. 

6. We have studied gas disk models which have 
density varying as a power law, truncated within 
an adjustable scale length, always preserving the 
Rayleigh stability of the disk, and a total mass 
from a fraction to dozens of Earth masses. The 
outcome of grain migration depends on the radia- 
tion pressure coefficient (3 and the grain size 5 . The 

5 The primary dependence is on (3. The approximate ranges 
of grain radii arc cited below for f3(s) based on our physical 
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fine dust fraction (/? > 0.5; s < 5/um) is removed 
by strong radiation pressure. Dust grains with 
0-1 < P < 0.5, or alternatively 5 < s < 20/zm, 
form an extended, outer dust disk component, re- 
siding most of the time beyond the bulk of the 
gas disk on eccentric orbits. Grains of intermedi- 
ate size (0.01 < < 0.1 or 20jan < s < 200(im) 
concentrate near the outer edge of the gas disk, 
forming a prominent dust ring of width similar to 
the scale length of gas density cutoff. Sand grains 
and meteoroids (0 < 0.01 or s > 200yum) stay in 
the inner part of the gas disk. 

8.2. Dust Morphologies and Their Ob- 
servability 

We considered collisions between each grain size 
fraction. The steady orbits can be called parking 
orbits at best. After one collisional time (which 
is typically much longer than the migration time, 
but much shorter than the age of the system) the 
grains are fragmented or eroded, jumping into a 
smaller grain fraction and moving to a new park- 
ing orbit or out into the interstellar space. The 
grains are continuously replenished, in our models 
from the inner part of the disk at radii of order 
10 AU, although a wider distribution of sources 
and/or somewhat different size distribution pro- 
duced by dust source are certainly thinkable. As 
a result, the disk reaches a stable dust density (and 
size) distribution. Including those particles which 
occupy the long-lived parking orbits, we predicted 
the dust distribution in several model disks. The 
results show an outer ring of small dust and, in ad- 
dition, an inner disk of large dust and sand. These 
results are somewhat sensitive to the mechanical 
properties of the particles, such as strength against 
catastrophic fragmentation and erosion coefficient. 
In contrast, the dust density in Figure 9 does not 
depend on the dust production rate in the disk. 

A variety of dust morphologies can be obtained 
in our models with deliberately smooth (mono- 
tonic) gas distributions, depending on how wide 
and dense the outer ring is, relative to the main 
disk. For instance, there might exist systems with 
ring morphology, ring plus inner disk morphology 

model of grains as moderately porous mixture of olivines 
and pyroxenes, and for the luminosity to mass ratio of the 
star equal 8.4, as appropriate for HR 4796A and approxi- 
mately also for HD 141569A. 



(possibly giving impression of disks with gaps), 
and disks with weak/unobservable rings. Since 
it is straitforward to obtain a ring several times 
denser than the adjacent disk region, such rings 
should be observable in scattered light imaging. 
Somewhat more model-dependent is the observ- 
ability of radial structure in thermal IR maps, 
since those depend strongly on the profile of dust 
temperature as a function of radius. 

We have qualitatively compared our models 
with two observed transitional disks, around HR 
4796A and HD 141569A (by transitional we mean 
object intermediate between the gas-rich primor- 
dial nebulae and late-stage, replenished dust disks 
of Pic and Vega-type; cf. the discussion of this 
notion in the introduction). We found a qualita- 
tive agreement with most of the observed charac- 
teristics of these disks, and sketched in §7.2 several 
possible improvements to the modeling procedure, 
which could be used in detailed model fitting. 

In general, our results open the prospect for di- 
agnosing the state of the gas component of circum- 
stellar disks by observing the dust and calculating 
the dust-gas coupling. The structure of the dust 
disk strongly depends on the properties of the gas 
disk. The most obvious example may be the ex- 
istence, radius and width of the outer "edge" or 
downturn of the gas disk, reflected in the corre- 
sponding parameters of the dust ring (of course, 
what we call a "ring" may well be in-out asymmet- 
ric, and thus prompt a description as a truncated 
disk). Therefore, if our model assumptions are 
satisfied, the dust component of an observed disk 
carries valuable information about the gas com- 
ponent, which is often more difficult or impossi- 
ble to study directly. For example, unlike in HD 
141569A, where the hydrogen and CO gas is at 
least directly detected if not mapped in detail, in 
the HR 4796A disk we only know the upper limits 
on gas contents. 

8.3. Origin of Structure in Disks: Dust 
Migration vs. Extrasolar Planets 

Our results also bear on the timely issue of 
planet detection in extrasolar systems. Radial and 
azimuthal structure can be induced in circumstel- 
lar disks by the gravitational perturbation (also 
known as sculpting) due to planets or other sub- 
stellar companions in or in the vicinity of disks 
(cf. Dermott et al. 1994, Liou & Zook 1999, Wy- 
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att et al. 1999, and Ozernoy et al. 2000, for 
gas- free disks and, e.g., Lubow et al. 1999 for 
gas disks). Detection of such structure is one of 
the main diagnostic tools in many of the future 
imaging projects at various wavelengths, includ- 
ing the ambitious space interferometry missions. 
However, caution should be exercised while inter- 
preting disk morphology. There are many possi- 
ble reasons unrelated to dynamical perturbations 
which could, in principle, masquerade as planets 
(Artymowicz 2000). 

First, optically thick disks with radius-dependent 
thickness might cast shadows on parts of them- 
selves. Next example would be a recent collision 
of planetesimals in disk. Though the cloud of 
dust released in such a collision would not be 
observable, it might become observable if expo- 
nentially amplified in the so-called dust avalanche 
(Artymowicz 1997, 2000). This requires the op- 
tical thickness in the disk mid-plane larger than 
that in the [3 Pic disk (transitional disks satisfy 
this requirement, as a rule). While admittedly 
rare and improbable, such a collision might pro- 
duce dense patches or sectors in a disk, which are 
non-axisymmctric and time-dependent on orbital 
timescale, features similar to those offered by dy- 
namical sculpting. Another example is the dust 
migration, potentially much more common, but 
fortunately easier to distinguish from planetary 
influences. 

Dust migration is capable of sculpting the 
disks, with an important difference from dynam- 
ical sculpting that it cannot itself produce non- 
axisymmctry (nor, presumably, any observable 
time variability). It can, however, provide an al- 
ternative explanation of some observed features, 
especially the radially structured dust disks, rings, 
and apparent gaps in disks. By carefully analyzing 
the axisymmetry (or lack thereof) of the under- 
lying dust distribution, taking properly into ac- 
count the asymmetry of the image due to the non- 
isotropic scattering function of a collection of dust 
grains, the observers will be able to decide whether 
or not our intrinsically axisymmetric models apply 
to any observed system. If axisymmetric model 
reproduce the observations, to within the uncer- 
tainties of the latter, and the inferred gas amount 
does not clash with any observational constraints, 
then the dust migration hypothesis is preferred 
over the planetary hypothesis by the virtue of be- 



ing more natural (conservative) . This is especially 
true in systems where gas has been independently 
detected. 

It appears that a convincing planet detection by 
way of dust imaging should include not just a non- 
axisymmetric (and time-dependent) structure in 
the disk, but a very specific one, closely following 
the predictions of numerical simulations of either 
gas-free or gas-rich disks. 
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Table 1 
Model Parameters. 



Model 


Star 


M g (M ffi ) 


rout(AU) 


V 


Cout 


K 


M d (M e ) 


L d /L 


1 


HR4796A 


10 


100 


-2.25 


1.05 


1 x 10~ 9 


0.90 


1.7 x 10~ 2 


2 




1 










0.10 


2.0 x 10~ 3 


3 








-2.75 






0.37 


8.2 x 10~ 3 


4 








-1.25 






3.10 


4.2 x 10~ 2 


5 










2.0 




0.72 


1.4 x 10" 2 


6 












2 x 10~ 8 


0.05 


9.1 x 10~ 4 


7 




4 


80 








0.31 


8.2 x 10~ 3 


8 


HD141569A 


50 


370 








4.23 


1.1 x 10" 2 



Note. — Blanks mean that the same value as model 1 is used. See §6.1 for the definition of K and §7.1 for 
M(j and L d . 
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Fig. 1. — Radiation pressure to gravity ratio (3 
around HR 4796A (solid line) and HD 141569A 
(dashed line). 
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Fig. 2. — Velocities of dust grains at 70 AU in 
the gas disk of model 1 described in §4. The cir- 
cles show the radial velocities v r and the diamonds 
show the relative velocities in the azimuthal direc- 
tion between the grains and gas, Ave — ve — v g . 
Filled symbols represent positive v r and negative 
Ave (outward migration) and open symbols repre- 
sent negative v r and positive Ave (inward migra- 
tion). The solid lines show the approximate for- 
mulae (23) and (25) calculated numerically after 
the iterations for the supersonic correction. The 
dashed lines show the first guesses of iteration, 
which are calculated using the subsonic stopping 
time T ss instead of T s . The horizontal dotted lines 
show the Kepler and thermal velocities at 70 AU 
in the model 1 disk. 
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Fig. 3. — Pressure gradient force to gravity ratio 
rj in the disk of model 1 (soled line) and model 5 
(dashed line). The disk in model 5 has wider edge 
at 100 AU than that in model 1. 
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Fig. 4. — Size of grains on the steady orbits with 
radii r. The solid and dashed lines show the grain 
size in model 1 and in model 5, respectively. The 
dotted line represents the case where the Poynting- 
Robertson (PR) drag is neglected in model 1 , and 
shows the difference from the solid line (including 
the PR drag) only for r > 130AU. 
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Fig. 5. — - (a) Evolution of the orbital radii r of 
8 — 300/im grains in the disk of model 1 . The mass 
of the gas disk is 10M ffi . The time is measured 
by the orbital period at 100AU, P ou t- The grains 
are generated at 10 AU and move outward until 
arriving at the steady orbits. The oscillations of r 
of 8[im and 20^m grains are due to the excitation 
of their orbital eccentricities at the edge of the gas 
disk, (b) The same as Fig. 5(a), but for model 2 
where the gas disk is !M ffi . 
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Fig. 6. — Evolution of the orbital elements of (a) 
an 8/im grain and (b) a 100/im grain. The gas disk 
is that in model 1. Arrows show the directions of 
the orbital evolution obtained by the perturbation 
equations (Gauss's equations). The length of ar- 
rows is set to be constant and does not represent 
the evolution rate. The circles show the evolu- 
tion of orbital elements of test particles, which 
are obtained directly by the orbital calculations. 
The initial values of the test particles' orbital el- 
ements are (a) a = 100AU and e = 0.05 and (b) 
a = 100AU and e = 0.5. The circles are plotted 
until t = lOOOPout with the logarithmic time in- 
terval: the first time interval is 0.01P ou t and each 
time interval is 1.2 times the previous one. 
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Fig. 7. — Distances of the apocenter and pericen- 
ter of the grains' orbits from the central star. The 
grains are generated at 10 AU and the orbits at 
t = lOOOPout are shown by the solid line (in the 
gas disk with 10M®; model 1) and dashed line 
(in the gas disk with 1M®; model 2). The orbits 
of grains larger than about 10/xm are circular and 
the distances of the apocenter and pericenter coin- 
cide with the value derived by the condition (3 = 7] 
(dotted line). 
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Fig. 8. — Life and migration time of dust grains. 
The solid line shows the life time of grains on 
the steady orbits under the collisional destruction. 
The dashed line shows the migration time from 
10AU to the 99% of the radius of the steady or- 
bits, (a) The gas mass is M g = 10M® (model 1). 
(b) M g = 1M® (model 2). 
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Fig. 9. — Optical depth in the vertical direction 
of the dust disks. The optical depth is calculated 
using the geometrical cross sections of grains, (a) 
For various disk masses, M g . the solid and dashed 
lines show the optical depth of the dust in the 
gas disk of 10M ffi (model 1) and lM e (model 2), 
respectively, (b) For various profiles of gas den- 
sity, p, where p g oc r p . The dotted, solid and 
dashed lines correspond to p = —1.25 (uniform 
surface gas density; model 4), p = —2.25 (model 
1) and p = —2.75 (model 3), respectively, (c) 
For various width of the edge of gas disks. The 
solid and dashed lines correspond to the narrow 
edge (Cout = 1.05; model 1) and the wide edge 
(Cout = 2.0; model 5), respectively, (d) For vari- 
ous strength of dust grains. The solid and dashed 
lines correspond to the rocky grains (K = 1 x 10 -9 ; 
model 1) and the icy grains (K = 2 x 10~ 8 ; model 
6), respectively. 
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Fig. 9 — continued. 
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Fig. 10. Simulated image of the model 4 disk 
in the scattered light and the thermal emission. 
The gray scale shows the scattered light flux in 
1.1/zm with linear spacing. The contours show the 
thermal flux in 18.2/zm and are spaced linearly. 
The unit of contour levels is arbitrary. The disk 
has a constant surface gas density with an edge 
at 100 AU, viewed in 60° from pole-on. The flux 
is calculated for r > 13AU (outside the artificial 
hole at the center). 
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